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Abstract. We show how estimates of parameters characterizing inflation- 
based theories of structure formation localized over the past year when 
large scale structure (LSS) information from galaxy and cluster surveys 
was combined with the rapidly developing cosmic microwave background 
(CMB) data, especially from the recent Boomerang and Maxima bal- 
loon experiments. All current CMB data plus a relatively weak prior 
probability on the Hubble constant, age and LSS points to little mean 
curvature (fltot = 1-08 ± 0.06) and nearly scale invariant initial fluctu- 
ations (n s = 1.03 ± 0.08), both predictions of (non-baroque) inflation 
theory. We emphasize the role that degeneracy among parameters in 
the L p k = 212 ± 7 position of the (first acoustic) peak plays in defin- 
ing the £ltot range upon marginalization over other variables. Though 
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the CDM density is in the expected range (O c d m h 2 = 0.17 ± 0.02), the 
baryon density f^h 2 = 0.030 ± 0.005 is somewhat above the independent 
0.019 ± 0.002 nucleosynthesis estimates. CMB+LSS gives independent 
evidence for dark energy = 0.66 ±0.06) at the same level as from su- 
pernova (SN1) observations, with a phenomenological quintessence equa- 
tion of state limited by SN1+CMB+LSS to w Q < -0.7 cf. the wq=-1 
cosmological constant case. 



1. CMB Analysis of Primary Anisotropies 

Experiments and Bandpowers: Anisotropies at the 30[iK level at low mul- 
tipoles revealed by COBE in 1992 were augmented at higher i in some 19 other 
experiments, some with a comparable number of resolution elements to the 600 
or so for COBE, most with many fewer. A list of these experiments to April 1999 
with associated bandpowers is given in Bond, Jaffe and Knox (2000 [BJK00]). 
The anisotropy picture dramatically improved this past year, as results were an- 
nounced first in summer 99 from the ground-based TOCO experiment in Chile 
(Miller et al. 2000), then in November 99 from Boomerang-NA, the North 
American test flight (Mauskopf et 1999). These two additions improved peak 
localization and gave evidence for Sltot ~ 1- Then in April 2000, results from the 
first CMB long duration balloon (LDB) flight, were announced (de Bernardis 
et al. 2000), followed in May 2000 by results from the night flight of Maxima 
(Hanany et al. 2000). Boomerang's best resolution was 10', about 40 times 
better than that of COBE, with tens of thousands of resolution elements. Max- 
ima had a similar resolution but covered an order of magnitude less sky. Fig. 1 
shows the 150A GHz Boomerang-LDB map and the Wiener-filtered Maxima-1, 
to scale. The de Bernardis et al. (2000) maps at 90 and 220 GHz show the 
same spatial features as this 150 GHz one, with the overall intensities falling 
precisely on the CMB blackbody curve. The Toco, Boomerang and Maxima 
experiments are described elsewhere in these proceedings. They were designed 
to reveal the primary anisotropies of the CMB, those which can be calculated 
using linear perturbation theory. Fig. 1 shows the temperature power spectra 
for Boomerang, Maxima and prior-CMB data (Boomerang-NA+TOCO+April 
99) are in good agreement. Sketching the impact of these new results on cosmic 
parameter estimation (Lange et al. 2000 [LetOO], Jaffe et al. 2000 [JetOO]) is 
the goal of this paper. Space constraints preclude adequate referencing here, 
but these are given in the Boomerang (LetOO) and Maxima+Boomerang (JetOO) 
parameter estimation papers (see also Bond 1996, [B96], for other references). 

We are only at the beginning of the high precision CMB era for primary 
anisotropies heralded by the arrival of Boomerang and Maxima, with interfer- 
ometers taking data (VSA, CBI, DASI), the single dish ACBAR about to, and 
new LDBs to fly in the next few years (Arkeops, Tophat, Beast/Boost), as well 
as Boomerang-2001 and the neo-Maxima Maxipol, both concentrating on polar- 
ization. In April 2001, NASA's HEMT-based MAP satellite will launch, with 
12' resolution, and in 2007, ESA's bolometer+HEMT-based Planck satellite is 
scheduled for launch, with 5' resolution. 
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Figure 1. The top figure shows Ci grouped in bandpowers for 
Boomerang-LDB (crosses), MaximaT (triangles) and prior-CMB experi- 
ments (TOCO+Boomerang-NA+" April 99", squares). The lower panel con- 
trasts the optimally-combined power spectra for Boomerang+Maxima+DMR 
(squares) with that for Boomerang+Maxima-|- prior-CMB (circles), showing 
the prior experiments do not move Ce very much. Best-fit models for arbi- 
trary Qtot and for f2 tot =l are shown in both panels. The Boomerang 150 A 
GHz map (i.e., for one of 16 bolometers) and the multifrequency Wiener- 
filtered Maxima-I map, its 124 square degrees drawn to scale, are shown in 
the bottom figure. Only the 440 square degrees within the central rectangle 
of the entire 1800 square degrees covered by Boomerang were used in the 
analvsis. 
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The CMB Analysis Pipeline: Analyzing Boomerang and other experiments 
involves a pipeline that takes (1) the timestream in each of the bolometer chan- 
nels coming from the balloon plus information on where it is pointing and turns 
it into (2) spatial maps for each frequency characterized by average tempera- 
ture fluctuation values in each pixel (Fig. 1) and a pixel-pixel correlation matrix 
characterizing the noise, from which various statistical quantities are derived, 
in particular (3) the temperature power spectrum as a function of multipole 
(Fig. 1), grouped into bands, and two band-band error matrices which together 
determine the full likelihood distribution of the bandpowers (Bond, Jaffe & Knox 
1998 [BJK98], BJKOO). Fundamental to the first step is the extraction of the sky 
signal from the noise, using the only information we have, the pointing matrix 
mapping a bit in time onto a pixel position on the sky. To compare the data with 
millions of cosmological models, as we wish to do here, the radical compression 
step from 2 to 3 is essential, and hinges upon an accurate representation of the 
likelihood surface. 

There is generally another step in between (2) and (3), namely separating 
the multifrequency spatial maps into the physical components on the sky: the 
primary CMB, the thermal and kinematic Sunyaev-Zeldovich effects, the dust, 
synchrotron and bremsstrahlung Galactic signals, the extragalactic radio and 
submillimetre sources. The strong agreement among the Boomerang maps in- 
dicates that to first order we can ignore this step, but it has to be taken into 
account as the precision increases. The Fig. 1 map is consistent with a Gaussian 
distribution, thus fully characterized by just the power spectrum. Higher or- 
der (concentration) statistics (3,4-point functions, etc.) tell us of non-Gaussian 
aspects, necessarily expected from the Galactic foreground and extragalactic 
source signals, but possible even in the early Universe fluctuations. For exam- 
ple, though non-Gaussianity occurs only in the more baroque inflation models 
of quantum noise, it is a necessary outcome of defect-driven models of structure 
formation. (Peaks compatible with Fig. 1 do not appear in non-baroque defect 
models, which now appear unlikely.) Though great strides have been made in 
the analysis of Boomerang and Maxima, there is intense needed effort worldwide 
now to develop new fast algorithms to deal with the looming megapixel datasets 
of LDBs and the satellites (e.g., Bond et al. 1999, Szapudi et al. 2000). 

2. Cosmic Parameters 

Parameters of Structure Formation: We usually adopt the restricted set of 
7 cosmological parameters used in LetOO and JetOO, ^b, ^cdm, n s , tc, os}- 

The curvature energy is f2fc = 1 — £ltot- The dark energy parameterized here by 
SI A could have complex dynamics associated with it, e.g., if it is the energy 
density of a scalar field which dominates at late times (now often termed a 
quintessence field, Q, with energy Qq, e.g., Steinhardt 2000). One popular phe- 
nomenology is to add one more parameter, wq = Pq/pq, where pq and pq are 
the pressure and density of the Q-field. Thus wq = — 1 and = £l\ for the 
cosmological constant. We have also allowed wq to float. 

We use 2 parameters to characterize the early universe primordial power 
spectrum of gravitational potential fluctuations <£, one giving the overall power 
spectrum amplitude V$(k n ), and one defining the shape, a spectral tilt n s (k n ) = 
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1 + d In V<& /din k, both at some (comoving) normalization wavenumber k n . We 
really need another 2, VGw(k n ) and rit(k n ), associated with the gravitational 
wave component. In inflation, the amplitude ratio is related to n t to lowest 
order, with 0(n s —n t ) corrections at higher order, e.g., B96. There are also useful 
limiting cases for the n s — n t relation. However, as one allows the baroqueness of 
the inflation models to increase, one can entertain a plethora of power spectra 
(with fully fc-dependent n s (k) and nt(k)) if one is artful enough in designing 
inflaton potential surfaces. As well, one can have more types of modes present, 
e.g., scalar isocurvature modes ("Pj s (A; n ), rii s {k)) in addition to, or in place of, 
the scalar curvature modes ("P$(fc n ), n s (k)). However, our philosophy is consider 
minimal models first, then see how progressive relaxation of the constraints 
on the inflation models, at the expense of increasing baroqueness, causes the 
parameter errors to open up. For example, with COBE-DMR and Boomerang, 
we can probe the GW contribution, but the data are not powerful enough to 
determine much. Planck can in principle probe the gravity wave contribution 
reasonably well. 

We use another 2 parameters to characterize the transport of the radiation 
through the era of photon decoupling, which is sensitive to the physical density 
of the various species of particles present then, Uj = Qjh 2 . We really need 4: ui, 
for the baryons, to c dm for the cold dark matter, LVhdm for the hot dark matter 
(massive but light neutrinos), and uj er for the relativistic particles present at that 
time (photons, very light neutrinos, and possibly weakly interacting products of 
late time particle decays). For simplicity, though, we restrict ourselves to the 
conventional 3 species of relativistic neutrinos plus photons, with uj er therefore 
fixed by the CMB temperature and the relationship between the neutrino and 
photon temperatures determined by the extra photon entropy accompanying 
e + e~ annihilation. Of particular importance for the pattern of the radiation is 
the (comoving) distance sound can have influenced by recombination (at redshift 
z r = a' 1 - 1), 

_ 6000 rV^ a\/a 

Ts ~ \/3 P °io (w m + ^ er a- 1 )V2(i + ^ a /( 4w7 /3)) 1 /2 ' i!; 

where c<j 7 = 2.46 x 10~ 5 is the photon density, uj er = 1.68a; 7 for 3 species of 
massless neutrinos and uj m = uJhdm + ^cdm + ^b- 
The angular diameter distance is 

TZ = {d k smh.(xr/d k ),Xr,d k sm(xr/d k )}, where d k = -^H== Mpc, (2) 



Xr = 6000 



M P C / t *r tttt 



The 3 cases are for negative, zero and positive mean curvature, d k is the cur- 
vature scale and Xr is the comoving distance to recombination. The location 
of the first acoustic peak, Lp k ~ 0.7467r7?./r s (e.g., Efstathiou and Bond 1999, 
hereafter EB99), depends upon uo\y through the sound speed as well as on uj k , ua 
and uj m . Thus Lp k defines a functional relationship among these parameters, a 
degeneracy (EB99) that would be exact except for the integrated Sachs- Wolfe 
effect, associated with the change of <3? with time if or £l k is nonzero. 
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Our 7th parameter is an astro-physical one, the Compton "optical depth" 
tc from a reionization redshift z re h to the present. It lowers Cg by exp(— 2tq) at 
the high £'s probed by Boomerang. For typical models of hierarchical structure 
formation, we expect tq ^ 0.2. It is partly degenerate with o% and cannot be 
determined at this precision by CMB data now. 

The LSS also depends upon our parameter set. Here we use a set of (rela- 
tively weak) constraints on In a\ from cluster abundance data and on T + (n s — 
l)/2 from galaxy clustering data (B96, LetOO). cr| is a bandpower for density 
fluctuations on a scale associated with rare clusters of galaxies, 8h _1 Mpc, which 
we often use in place of V$>(k n ) or Cio for the amplitude parameter. The mass- 
density-power-spectrum-shape-parameter r depends upon {u> m , uj er , uj^, h} and 
is related to the horizon scale when the energy density in relativistic particles 
equals that in nonrelativistic ones. 

When we allow for freedom in uj er , the abundance of primordial helium, tilts 
of tilts (dn^ s is t y(k n )/dliik, ...) for 3 types of perturbations, the parameter count 
would be 17, and many more if we open up full theoretical freedom in spectral 
shapes. However, as we shall see, as of now only 3 or 4 combinations can be 
determined with 10% accuracy with the CMB. Thus choosing 7 is adequate 
for the present, 6 of which are discretely sampled, with generous boundaries. [J] 
For drawing cosmological conclusions we adopt a weak prior probability on the 
Hubble parameter and age: we restrict h to lie in the 0.45 to 0.9 range, and the 
age to be above 10 Gyr. 

The First Peak, Qtot and Lpk and its errors are found from the average 
and variance of In Lp^, taken wrt the full probability function over our database 
described above (restricted for this exercise to the tq = part and with the weak 
prior). As more CMB data were added, Lpk evolved and the errors shrunk con- 
siderably: 240i^ for April 99 data, 220±^ for TOCO+4.99 data, 224±|f when 
Boomerang-NA was added, and 212_ 7 when Boomerang-LDB and Maxima-1 
were added to prior-CMB. The latter contrasts with 202^£ for Boomerang-LDB 

alone, 226^Jg for Maxima-1 alone, and 208^y for the combination. The numbers 
change a bit depending upon exactly what prior one chooses or what functional 
forms one averages over. For the choice here, although there is a large difference 
in the mean between the Maxima and Boomerang numbers, it is not unreason- 
able within the errors. Other ways of doing this make the discrepancy seem 
more statistically significant (e.g., Page 2000, these proceedings). 

In Fig. 2, we show the lines of constant Lp^ cx 1Z/r s in the Qtot~^A plane, 
for given ui m and u>b, using the formulas given above and discussed in more detail 
in EB99. The ±10 band around 210 corresponds to our best Lpk estimate of 



The specific discrete parameter values used for the C^-database in this analysis were: (Qa = 
0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0,1.1), (Q k = .9,.7,.5,.3,.2,.15,.l,.05,0,-.05,-.l,-.15,-.2,-.3,-.5), (r c =0, 
.025, .05, .075, .1, .15, .2, .3, .5) when wq = -1, with slightly different ranges (and fl tot — l) 
when we allow wq to float; (lu c = .03, .06, .12, .17, .22, .27, .33, .40, .55, .8), (w b = .003125, 
.00625, .0125, .0175, .020, .025, .030, .035, .04, .05, .075, .10, .15, .2), (n s =1.5, 1.45, 1.4, 
1.35, 1.3, 1.25, 1.2, 1.175, 1.15, 1.125, 1.1, 1.075, 1.05, 1.025, 1.0, .975, .95, .925, .9, .875, .85, 
.825, .8, .775, .75, .725, .7, .65, .6, .55, .5), erf was continuous, and there were 4 experimental 
parameters for Boomerang and Maxima (calibration and beam uncertainties), as well as other 
calibration parameters for some of the prior-CMB experiments. 
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212 ± 7 using all current CMB data. Note that the constant Lp^ lines look 
rather similar to the contours shown in the right panel, showing that the lZ/r s 
degeneracy plays a large role in determining the contours. The contours hug the 
^tot = 1 line more closely than the allowed Lp^ band does for the maximum 
probability values of uj m and u>b, because of the shift in the allowed Lpj~ band 
as uj m and lu^ vary in this plane. See also Bond et al. (2001b, [capp2K]) for Lpj, 
plots in the quintessence plane, wq-Qq, which demonstrate why wq is poorly 
determined by CMB alone. 



Table 1. Cosmological parameter values and their 1-sigma errors are 
shown, determined after marginalizing over the other 6 cosmological and 4 + 
experimental parameters, for B98+Maxima-I+prior-CMB and the weak prior, 
0.45 < h < 0.9, age > 10 Gyr. The LSS prior was also designed to be weak. 
In the first set fi tot varies, in the second set it is fixed to unity. Similar tables 
for B98+DMR are given in LetOO and for B98+MAXIMA-I+DMR in JetOO. 
We have set the quintessence wq = Pq/pq parameter to —1, the cosmolog- 
ical constant case, but the last line shows the limit on wq if we allow it to 
vary (the other parameters do not move much). SN1 results on wq with 
earlier CMB data were given in Perlmutter et al. (1999). The detections 
in the table are clearly very stable if extra "prior" probabilities for LSS and 
SN1 are included, and are also stable with much stronger priors on h, but 
do move if the BBN-derived 0.019 ± 0.002 prior is imposed. If Vt to t is varied, 
parameters derived from our basic 7 come out to be: age=13.2 ± 1.3 Gyr, 
h = 0.70 ± 0.09, VL m = 0.35 ± .06, n b = 0.065 ± .019. Restriction to fi tot = 1 
yields: age=11.6±0.4 Gyr, h = 0.80±.04, VL m = 0.31±.03, n b = 0.046±.005. 



-LSS 



-SN1+LSS 



unb 



-SN1 



TU7 - 

.07 
.005 



h.06 
-.06 
n oi+.005 
• u<3i -.005 

• i4 -.02 

1 04+- 09 
0-63±;g| 



.05 
.005 



-.04 
0T1+- 005 

i r; + -03 
• io -.02 

1.06l;gj 
0.70_^ 



Qtot 

h 



n h h 2 



n s 

n A 



1.09; 



•0311.005 

1 7+.O6 
— .05 

i.uo_ 08 
0-48±;|° 



1.08] 



1.04; 



1.04] 



•0311.005 

5+.05 
.05 
h.10 
-.09 
-.07 
-.07 



.13] 

Los: 

0.72 H 



n 



tot 



--1 



CASE 



DOT" 
004 
1O+-06 
• iy -.05 

L.02±;gf 



■ U6U ~.004 
i 7+.02 
• i '-.02 

1.03tf 7 



.uau_ 004 
.16l° 3 
1.03 



.UJU_ 004 

17 +.oi 

■ L '-.02 
.07 
-.07 
.03 
.05 



n b h 2 

n s 

n A 



mo: 



0.711 



wq (95%) 



.03 
f .08 
-.07 
-.06 
.07 



1.04 
0.69 



< -0.29 < -0.33 < -0.69 < -0.73 



Marginalized Estimates of our Basic 7 Parameters: Table 1 shows there 
are strong detections with only the CMB data for fttot, w 6 an d n s in the min- 
imal inflation-based 7 parameter set, and a reasonable detection of Lo c dm- 

The 

ranges quoted are Bayesian 50% values and the errors are 1-sigma, obtained after 
projecting (marginalizing) over all other parameters. That Cl\ is not well deter- 
mined is a manifestation of the Qtot~^A near-degeneracy discussed above, which 
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dark matter + baryons dark matter + baryons 




L Pk for cj m =.17,w b =.03,w (! =-l L Pk for w m = .25,cj b =.02,w Q =- 1 




Figure 2. The top panels show 1,2,3-sigma likelihood contours for the 
weak-H+age prior probability (right) and when the LSS "prior" is included 
(left). In the right panel, "prior-CMB" experiments (TOCO+Boomerang- 
NA+" April 99") were included instead of just DMR, but the figure is very 
similar for Boomerang+Maxima+DMR. The supernova contours are also plot- 
ted, and the solid contour lines are what you get when you combine the two 
likelihoods. The bottom panels show lines of constant Lpk in the £l m —£lA 
plane for two choices of {w m ,W(,}, left the most probable values, right when 
the current BBN constraint is imposed (lowering increases the sound speed, 
decreasing Lpk, and varying uj m also shifts it). The 0.65 < h < 0.75 (heavier 
shading, H ) and 11 < age < 15 (lighter shading, to) ranges and decelerations 
<7o = 0, —1/3, —1/2 are also noted. The sweeping back of the Lpk curves into 
the closed models as ft a is lowered shows that even if J7t ot =l is the correct 
answer, the phase space results in a ID projection onto the f2 tot axis that 
would be skewed to f2 tot > 1, a situation we see in Table 1. Note that the 
contours in the top left panel are near the diagonal £ltot = 1 line, but also 
follow a weighted average of Lpk ~ 210 lines. This approximate degeneracy 
implies 57a is poorly constrained for CMB-only, but it is broken when LSS 
is added, giving a solid SNl-independent Qa "detection". Lpk and contour 
plots in the wq-Qq plane for O tot =l are given in capp2K. 
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is broken when LSS is added because the CMB-normalized as is quite different 
for open cf. A models. Supernova at high redshift give complementary informa- 
tion to the CMB, but with CMB+LSS (and the inflation-based paradigm) we 
do not need it: the CMB+SN1 and CMB+LSS numbers are quite compatible. 
In our space, the Hubble parameter, h = -(fijh 2 )) 1 / 2 , and the age of the 

Universe, to, are derived functions of the Qjh 2 : representative values are given 
in the Table caption. 

Fig. 3 shows how the parameter estimations evolved as more CMB data 
were added (for the weak+LSS prior). With just the COBE-DMR+LSS data, 
the 2-sigma contours were already localized in oo c dm- Without LSS, it took the 
addition of Maxima-1 before it began to localize. £lk localized near zero when 
TOCO was added to the April 99 data, more so when Boomerang-NA was added, 
and much more so when Boomerang-LDB and Maxima-1 were added. Some n s 
localization occurred with just "prior-CMB" data, uj^ really focussed in with 
Boomerang-LDB and Maxima-1, as did 

We have also considered what happens as we let Q, mi/ /Q, m , the fraction of the 
matter in massive neutrinos, vary from to 0.3, for LSS + all of the CMB data 
and n tot =l (Bond et al. 2001a, [v2K]). Until Planck precision, the CMB data 
by itself will not be able to strongly discriminate this ratio. Adding HDM does 
have a strong impact on the CMB-normalized as and the shape of the density 
power spectrum (effective T parameter), both of which mean that when LSS 
is included, adding some HDM to CDM is strongly preferred in the absence of 
JIa- However, though higher f2 m is preferred at the expense of less dark energy, 
significant f^A is still required (see v2K for the evolution of the CMB+LSS 2- 
sigma contours in the {oJhdm + ^>cdm)^A plane as Q. mv /Q. m is varied). The cob 
and n s likelihood curves are essentially independent of Q mu /£l m . 

The Future, Forecasts for Parameter Eigenmodes: We can also forecast 
dramatically improved precision with further analysis of Boomerang and Max- 
ima, future LDBs, MAP and Planck. Because there are correlations among the 
physical variables we wish to determine, including a number of near-degeneracies 
beyond that for Q tot -0,\ (EB99), it is useful to disentangle them, by making 
combinations which diagonalize the error correlation matrix, "parameter eigen- 
modes" (e.g., B96, EB99). For this exercise, we will add uJhdm and n t to our 
parameter mix, but set wq=— 1, making 9. (The ratio VGw(k n ) /Vq>(k n ) is 
treated as fixed by n^, a reasonably accurate inflation theory result.) The fore- 
cast for Boomerang+DMR based on the 440 square degree patch with a single 
150 GHz bolometer used in the published data is 3 out of 9 linear combina- 
tions should be determined to ±0.1 accuracy. This is indeed what we get in 
the full analysis of LetOO for CMB only. If 4 of the 6 150 GHz channels are 
used and the region is doubled in size, we predict 4/9 could be determined to 
±0.1 accuracy. The Boomerang team is still working on the data to realize this 
promise. And if the optimistic case for all the proposed LDBs is assumed, 6/9 
parameter combinations could be determined to ±0.1 accuracy, 2/9 to ±0.01 
accuracy. The situation improves for the satellite experiments: for MAP, we 
forecast 6/9 combos to ±0.1 accuracy, 3/9 to ±0.01 accuracy; for Planck, 7/9 to 
±0.1 accuracy, 5/9 to ±0.01 accuracy. While we can expect systematic errors to 
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Figure 3. 2-sigma likelihood contours for the dark matter density lo c — 
^cdmh 2 arid {fife, f^A, n s , Wf,} are plotted for LSS, the weak-H+age cosmo- 
logical prior, and the following CMB experimental combinations: DMR 
(short-dash); the "April 99" data (short-dash long-dash); TOCO+4.99 data 
(dot short-dash); Boomerang-NA+TOCO+4.99 data (dot long-dash, termed 
"prior-CMB"); Boomerang-LDB+Maxima-l+Boomerang-NA+TOCO-|-4.99 
data (heavy solid, all-CMB). These 2a lines tend to go from outside to inside 
as more CMB experiments are added. The smallest 2-sigma region (dotted 
and interior) shows SNl+LSS+all-CMB, when SNI data is added. For the 
f^A, n s and uj^ plots , Q to t—^ has been assumed, but the values do not change 
that much if Q, to t floats, as Table 1 shows for all-CMB. (The discreteness of 
the C^-database is responsible for the sharpness of the contour edges, espe- 
cially evident in the lo c -VI^ CMB+LSS case. Though different interpolation 
schemes can round these off somewhat, it is the price of our finite grid.) 
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loom as the real arbiter of accuracy, the clear forecast is for a very rosy decade 
of high precision CMB cosmology that we are now fully into. 
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